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ABSTRACT 

Observations of the central radian of the Galaxy by the Reuven Ramaty 
High Energy Solar Spectroscopy Imager (RHESSI) have yielded a high-resolution 
measurement of the 1809 keV line from 26 Al, detected at llcr significance in 
nine months of data. The RHESSI result for the width of the cosmic line is 
(2.03±5;|?) keV FWHM. The best fit line width of 5.4 keV FWHM reported by 
Naya et al. (1996) using the Gamma-Ray Imaging Spectrometer (GRIS) balloon 
instrument is rejected with high confidence. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — line:profiles 
— gamma rays:observations 



1. Introduction 

The 1809 keV line of 26 A1 was the first astronomical gamma-ray line from radioactive 
decay discovered in the Galaxy (Mahoney et al. 1984), confirming a prediction made by 
Ramaty & Lingenfelter (1977) and Arnett (1977). It is the brightest Galactic line after the 
line from positron annihilation at 511 keV. This first result was from HE AO 3, which used 
high-resolution germanium detectors. Mahoney et al. (1984) found that the line width was 
consistent with their instrumental resolution of about 3 keV, but the statistical significance 
of the result was low (4.8 a). 

This isotope is thought to be produced by proton capture on magnesium, and therefore 
can occur in any site of nucleosynthesis where these two components are abundant (Prantzos 
& Diehl 1996). Environments suggested for its creation and release into the interstellar 
medium include type-II supernovae, novae, and the winds of Wolf-Rayet and Asymptotic 
Giant Branch stars. Because its half-life is around 10 6 yr, much shorter than scales of 
Galactic evolution, its distribution on the sky reflects the current Galactic distribution of 
the relevant parent events. Maps of the Galactic 1809 keV emission (Knodlseder et al. 
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1999, and references therein) were made with the COMPTEL instrument on the Compton 
Gamma-Ray Observatory and have given us the first detailed look at where the emission is 
concentrated. The maps correlate better with early than late stellar populations, suggesting 
that supernovae and/or Wolf-Rayet stars are likely to be the primary contributors. 

Many other observations of this line, with high- and low-resolution balloon instruments 
and with low-resolution satellite instruments, have occurred since HE AO 3 (see Prantzos & 
Diehl (1996) for a comprehensive review). Only one high-resolution measurement has had a 
statistical significance comparable to the HE AO 3 data and thus been able to further advance 
our understanding of the shape of the line. This was performed by the Gamma-Ray Imaging 
Spectrometer (GRIS) balloon (Naya et al. 1996). The authors found the line to be signifi- 
cantly broadened, with an intrinsic width of (5. 41^) keV derived from a measured width of 
(6.4+J; 2 ) keV by subtracting their instrumental width of (3.4 ± 0.1) keV in quadrature. The 
significance of their overall detection was 6.8a, slightly higher than that of HEAO 3. 

This Doppler broadening corresponds to isotropic velocities of 540 km s _1 or a temper- 
ature of ~ 4.5 x 10 8 K (Sturner & Naya 1999), and there is no model for how either could 
be maintained by a gas in the interstellar medium for anything approaching the lifetime of 
the isotope. This result has stimulated interesting theoretical work centered on concentrat- 
ing the 26 Al in grains, which can maintain their birth velocities much longer than gaseous 
material and even be re-accelerated in supernova shocks (Chen et al. 1997; Ellison, Drury, 
& Meyer 1997; Lingenfelter, Ramaty & Kozlovsky 1998; Sturner & Naya 1999). 

2. The instrument and analysis technique 

The Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) is a NASA 
Small Explorer satellite in a nearly circular ~600 km orbit (96-minute period) with inclina- 
tion ~38°. Its primary mission is to make high-resolution images and spectra of solar flares 
in the range from 3 keV to 17 MeV (Lin et al. 2002, 2003). Its detector array is a set of 
nine high-purity coaxial germanium detectors cooled to liquid nitrogen temperature by a 
Stirling-cycle refrigerator (Smith et al. 2002). Each of RHESSFs detectors is segmented into 
a thin front segment facing the Sun (to stop hard x-rays) and a much thicker rear segment 
meant to record solar gamma-rays that penetrate the front segment. Because the front seg- 
ments have little effective area at 1809 keV for any incident angle, I did not use them for 
this analysis. One of the nine detectors operates in an unsegmented mode with poor energy 
resolution, and is excluded as well. The spacecraft rotates at close to 15 rpm about the axis 
pointing at the Sun. 
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Because the array is unshielded and the spacecraft is very light, the effective area of 
the array to highly penetrating 1809 keV photons is nearly independent of the direction of 
incidence, either in azimuth or zenith angle with respect to the spacecraft axis. This has been 
verified with Monte Carlo simulations using the GEANT3 package with a highly detailed 
model of the spacecraft, which show a maximum deviation of ±10% at any angle from the 
mean effective area averaged over all angles. Thus the annual rotation of the spacecraft 
with respect to the stars cannot be used to modulate the Galactic signal as was done with 
data from the Solar Maximum Mission Gamma- Ray Spectrometer (Harris et al. 1990). As 
an additional complication, there is a background line in the instrument at 1809 keV due 
to cosmic-ray interactions with aluminum in the spacecraft. This must be understood and 
subtracted before the Galactic line can be studied. 

The lack of directionality and strong background line can be overcome by using the 
Earth as an occulter. To begin the analysis, I divided nine months of RHESSI data (1 
March 2002 to 2 December 2002) into one-minute (and therefore spin-averaged) intervals. I 
defined the "inner Galaxy" as a box running from ±30° in Galactic longitude and ±5° in 
Galactic latitude. "Source" intervals were times when this entire box was unocculted by the 
Earth, and "background" intervals were defined as times when it was entirely occulted. The 
rest of the data were discarded. Because the Earth subtends less than half the sky, the set 
of source pointings is larger. The total amount of time accumulated in the source spectra 
is 75.0 dy, and the total accumulation of background data is 33.7 dy. Data contaminated 
by the precipitation of magnetospheric electrons were identified with the onboard particle 
detector and were removed from consideration, as were data taken during the X4.8 solar flare 
of 23 July, the only flare observed so far with emission at energies approaching 1809 keV. 

The background subtraction algorithm is adapted from one used to derive the most 
precise value to date for the flux of positron-annihilation radiation from the Galactic Center 
region: (1.49 ± 0.02) xl0~ 3 ph cm~ 2 s _1 (Smith et al. 1997). In that case the instrument 
was the Burst and Transient Source Experiment (BATSE) on the Compton Gamma-Ray 
Observatory ( CGRO). The principle is to sort the background spectra into a two-dimensional 
library based on the two most important parameters that control the background. For each 
source spectrum, a background spectrum is generated by selecting among the library spectra 
using the values of the controlling parameters for that source spectrum. The computer code 
tracks the repeated contributions of each original raw background spectrum in order to 
propagate statistical errors correctly. 

All background lines with half-lives less than a few hours are fairly well controlled 
by choosing two particular parameters on which to sort the library: the longitude of the 
ascending node (LAN) of the current orbit, and the orbital phase since the ascending node 
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(PAN). These parameters define a particular spot on the Earth, so that the cosmic-ray flux 
(which is controlled primarily by geomagnetic latitude) is reproduced. The lines with half- 
lives of seconds or less are proportional to the instantaneous cosmic-ray flux: this includes the 
1809 keV background line due to cosmic-ray production of 26 *Mg in aluminum, which decays 
with a 0.49 ps half-life (Wheaton et al. 1989). Unlike simple geographic or geomagnetic 
coordinates, however, LAN and PAN also contain a memory of the recent history of the 
spacecraft trajectory - most usefully, how recently and how deeply the spacecraft passed 
through the SAA. Thus matching LAN and PAN also produces a good match for the SAA- 
induced background lines, which have half-lives from minutes to hours. We see such a line 
around 1811 keV that blends with the 1809 keV line. The only candidate in the gamma-ray 
tables of Chu, Ekstrom & Firestone (1999) is the 1810.77 keV line of 56 Mn (half-life 2.6 hr), 
presumably created by interactions of SAA protons with iron near the detectors. This line 
was suggested as a contributor to the background of the germanium spectrometer on the 
Wind spacecraft (Weidenspointner et al. 2001), and the published 1809 keV background line 
from HE AO 3 (Mahoney et al. 1984) shows an excess in the blue wing that is consistent with 
the 56 Mn line appearing at a low level. There is no very long-lived background line (such 
as the 1275 keV line from 22 Na) near 1809 keV. Such a line would require a background- 
selection system which is cognizant of the total time since the start of the mission, rather 
than relying entirely on LAN and PAN. 



3. Results 

Figure 1 (left) shows the RHESSI background-subtracted spectrum of the Galactic 26 Al 
line. The smooth curve shown is 3% of the average background for comparison. Fitting 
the spectrum from 1790-1825 keV to a Gaussian plus a linear background, I find a center 
energy of (1808.87±0.18) keV, 1.2<r from the expected rest energy of the line at 1808.65 keV 
(Chu, Ekstrom & Firestone 1999). The of the line is (1.17 ± 0.11) xlO" 2 counts s" 1 
(11<t), which is 12.9% of the background line (averaged over the entire set of background 
spectra interpolated from the background library). Its FWHM, which includes instrumental 
broadening, is (4.58 ± 0.44) keV. The fit is good (x 2 = 26.7 with 30 degrees of freedom). 

By fitting seven narrow background lines ranging from 186 keV to 2243 keV, I interpolate 
the intrinsic instrumental resolution to its value at 1809 keV and find (4.10 ± 0.07) keV. I 
avoid the instrumental background line at 1809 keV itself because of the blend with the 
line at 1811 keV, which makes it much broader than the interpolation ((4.86 ± 0.08) keV 
FWHM) and raises its center energy to (1809.34 ± 0.03) keV. Subtracting the instrumental 
resolution in quadrature from the Galactic line gives an intrinsic width of (2.03TJ?;2i) keV, 
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corresponding to velocities of ~200 km/s or a temperature of 1.7 x 10 8 K. 

To demonstrate the quality of the background subtraction, I show in Figure 1 (right) the 
region around a background line at 1014 keV. This line, like the background line at 1809 keV, 
is due to the prompt decay of a short-lived isomer created by cosmic-ray interactions in 
aluminum. Thus, if it subtracts well, one should be able to assume that the 1809 keV 
background line does also. The figure demonstrates subtraction better than the 3% shown 
for comparison. There is a hint of a residual 1014 keV line which, when fit with a Gaussian, 
gives a flux of (1.2 ± 0.5)% of the background line. The averaged cosmic-ray count rate in 
the detectors over all the periods used for source pointings is 0.46% higher than the same 
average over the background periods. This is consistent with the slight undersubtraction 
hinted at in Figure 1 (right). I find that subtracting out an extra 0.46% or 1.2% of the 
1809 keV background line has no significant effect on the Galactic 1809 keV line width, 
whether I use the natural form of the background spectrum (which includes the 1811 keV 
line) or whether I substitute an artificial single Gaussian line at exactly 1808.65 keV with 
the instrumental resolution. Since the Galactic line flux is 12.9% of the average background, 
these extra subtractions would reduce its value by 4% and 9%, respectively. 

To further demonstrate that the cosmic-ray-induced background line is well-subtracted, 
I divided the data set into two roughly equal parts by the value of the cosmic-ray count 
rate during each one-minute spectrum. The rate averaged 158 counts s" 1 for the lower data 
set and 237 counts s -1 for the higher data set, for a ratio of 1.50. The count rate in the 
total, unsubtracted 1809 keV line had a similar ratio of 1.42, with 8.00 xlO -2 counts s" 1 and 
11.36 xl0~ 2 counts s _1 in the two sets, respectively. Once the backgrounds were subtracted, 
however, the values for the residual (Galactic) line were (1.13 ± 0.15) xl0~ 2 counts s _1 and 
(1.23 ± 0.15) xl0~ 2 counts s -1 , respectively, a ratio of only 1.09 and well within statistical 
agreement. The derived Galactic signal thus appears to be independent of the intensity of 
the background line. 

The conversion of the count rate in the Galactic line to a total Galactic flux is strongly 
dependent on the assumed Galactic distribution. The effective area of the eight rear segments 
used, calculated at 1809 keV using the distribution of angles of incidence over the observation 
period assuming the source is concentrated at the Galactic center, is 20.5 cm 2 . The incident 
flux is then (5.71 ± 0.54) xl0~ 4 ph cm -2 s _1 for the artificial case of a point source. Since 
the effective area of the instrument is nearly independent of angle, the model distribution 
is important in determining the flux only because of occultation of parts of the distribution 
by the Earth. As noted above, I required the inner 60° of Galactic longitude to be visible in 
each pointing, but on average a considerable additional amount of the plane is also included. 
A certain amount of this high-longitude plane emission also appears in the background 
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intervals, however, and is subtracted off. 

4. Discussion 

Adding the GRIS best-fit value of 5.4 keV for the Galactic line width in quadrature 
with RHESSTs instrumental resolution would give a 6.78 keV width. Fixing the width at 
that value for the fit from 1790-1825 keV causes x 2 to increase by 17.3 to 44.0 (now with 
31 degrees of freedom). The probability of RHESSrs result being consistent with the GRIS 
best fit is then 4xl0~ 5 for one parameter of interest. 

The RHESSI value for the Galactic line's width is only marginally consistent with zero. 
It is expected, however, that the interstellar 26 Al will share in Galactic rotation and display 
appropriate Doppler shifts (Skibo & Ramaty 1991; Gehrels & Chen 1996). Gehrels & Chen 
(1996) combined a three-dimensional model 26 Al distribution derived from COMPTEL data 
with a Galactic rotation model to produce a map of the radial velocity versus Galactic 
longitude (integrated over Galactic latitude and line of sight). This map has a Galactic 
bulge component with velocities up to ~75 km/s and spiral arm components with velocities 
up to ~150 km/s. Although an integration of this map over the inner Galaxy has not been 
published, from the slices shown at different longitudes it appears that a integrated FWHM 
of not more than 150 km/s or about 0.9 keV would be obtained. This is within la of our 
result. Additional broadening due to more local motions of the 26 Al-bearing gas or dust is 
possible, and of course the most recently-created portion of the 26 Al can be traveling near 
its birth velocity in supernovae, etc., even if this portion is small. 

RHESSTs flux value of (5.71 ± 0.54) xl0~ 4 ph cm~ 2 s _1 is comparable to previous 
measurements. A value of ~ 4xl0~ 4 ph cur 2 s~ x rad -1 was considered consistent with 
the existing ensemble of data sets at the time by Diehl & Timmes (1997). For a uniform 
distribution in the plane, this would imply about 1.4 rad for RHESSFs effective field of view 
(including the effect of subtracted flux in the background pointings). Realistic distributions 
fall off with longitude, however, and would require an even larger effective field of view. 
Naya et al. (1998), however, using the GRIS data, found higher fluxes per radian, with 
the value depending on the model used and ranging up to (5.48 ± 0.78) xl0~ 4 ph cm~ 2 s _1 
rad -1 assuming the COMPTEL distribution available at that time. They suggested that 
COMPTEL, which subtracts background by imaging with the Compton telescope technique, 
would be insensitive to a spatially diffuse component that the large, integrating field of view 
of GRIS would accept. At first glance the RHESSI result seems to support this higher value, 
since the effective field of view required to match the GRIS data, 1.04 radians, is identical 
to the region I required to be visible for the source intervals. In future work, these RHESSI 
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source and background fields of view, with the addition of data from those times in which the 
central radian is partially occulted, will be convolved with different model distributions to 
provide a normalization for each distribution and to help choose among them. Comparison 
with upcoming results from INTEGRAL/SPI, which has a smaller field of view, will allow 
us to test for the existence of very large-scale diffuse emission. 

I would like to thank the rest of the RHESSI team for making this work possible, 
particularly those who worked on the spectrometer and the data- analysis software. G. H. 
Share made valuable suggestions related to this analysis, and R. P. Lin and H. Hudson 
thoughtfully reviewed the manuscript. This work was supported by NASA contract NAS5- 
98033. 
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Fig. 1. — Inner-Galaxy spectra in the vicinity of the 1809 keV line (left) and the background 
line at 1014 keV (right). In each plot, the lower curve is the RHESSI background-subtracted 
count spectrum and the upper curve is 3% of the average background spectrum during this 
period. The bins shown are 1 keV wide and the instrumental resolution at 1809 keV is 
(4.10 ±0.07) keV (see text). 



